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ABSTRACT 

We report variability of the X-ray source, X-7, in NGC 6946, during a 60 ksec Chandra observation 
when the count rate decreased by a factor of ~ 1.5 in ~ 5000 sees. Spectral fitting of the high and low 
count rate segments of the light curve reveal that the simplest and most probable interpretation is that 
the X-ray spectra are due to disk black body emission with an absorbing hydrogen column density 
equal to the Galactic value of 2.1 x 10^^ cm^^. During the variation, the inner disk temperature 
decreased from ~ 0.29 to ~ 0.26 keV while the inner disk radius remained constant at ^ 6 x 10^ cm. 
This translates into a luminosity variation from 3.8 to 2.8 x 10'^^ ergs cm~^sec~^ and a black hole 
mass of ~ 400Mq. More complicated models like assuming intrinsic absorption and/or the addition of 
a power-law component imply a higher luminosity and a larger black hole mass. Even if the emission 
is beamed by a factor of ~ 5, the size of the emitting region would be > 2.7 x 10^ cm implying a black 
hole mass > 18OM0. Thus, these spectral results provide strong evidence that the mass of the black 
hole in this source is definitely > 100 Af© and more probably ^ 4OOM0. 
Subject headings: accretion, accretion disks — galaxies: individual (NGC6946) — X-rays: binaries 



1. INTRODUCTION 

Chandra observations of nearby galaxies, have 
confirme d the preseii c e of n on-nuclear X-ray point 
sources ([K aarct et aT 12001 1 : iMatsumoto et al\ 120011: 
IZezas fc Fa bbiano 2002), which have luminosities > 10'^^ 
ergs/s and hence have been called Ultra luminous X-ray 
Sources (ULX). 

Since these sources radiate at a rate greater than 
the Eddington luminosity for a ten-solar mass black 
hole, they are believed to harbor a b lack hole of mass 
10 M© < M < lO^M © (|Colbert fc Mushotzkvl 119991 : 
iMakishima et al.] l2000( ) where the upper limit is con- 
strained by the argument that a more massive black hole 
would have s ettled into the nucleus due to dynamical 
friction (Kaa ret et al\ 12001 f ). If this interpretation is 
true, than these black holes have mass in the interme- 
diate mass range between those of stellar mass black 
holes found in Galactic X-ray binaries and those asso- 
ciated with Active Galactic Nuclei and hence are called 
Intermediate Mass Black Hole s (IM BH) . F o r a re view see 
iMiUer. Fabian fc Milled (pOOi ) andlMillerl pOOSh . 

The creation of s uch black holes 

(iPortegies. Zwart fc McMillianI l2002t iTaniguchi et al\ 
I2000t iMadau fc Reed 1200 It) and the proces s by which 
they sustain such high accretion rates (jKing et al\ 
|2001| ). are largely unknown and and when understood 
are expected to make radically shifts in the present 
paradigms of stellar and binary evolution and the history 
of the Universe. On the other hand, alternate models for 
ULX challenge our present understandi ng of accreting 
systems such as super-Eddington disks ( Begelmanl l2002l) 
or emission tha t is beamed from a geometrically thick 
accretion disk ([King et a/.ll2001t ). For the latter case, 
it has been argued that such thick "funnel" shaped 
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disks enhance the observed flux by just a factor of 
few (jMisra fc Sriram 2003). Thus, it is important to 
ascertain whether ULX do indeed harbor IMBH or not. 

Since a more direct measure of the mass such as spec- 
troscopic mass function measurement of the binary, is 
not possible for ULX, indirect evidences have to be used. 
One such way is to look for similarities in the spectral 
and temporal properties of ULX and black hole X-ray 
binaries. There are theoretical indications that the na- 
ture of the accretion flow should depend on the Edding- 
ton ratio L/Ledd rather than on the actual values of the 
bolometric luminosity L and the Eddington limit LEdd- 
Thus, a ULX should display similar spectral and tempo- 
ral characteristics as a black hole X-ray binary accreting 
at a similar L/ LEdd, even though the masses of the black 
holes are different. 

If ULX harbor IMBH, they should display analogues 
spec tral states to X-ra y binaries. Based on such a anal- 
ogy, lYuan et ahl (|2007D modeled the X-ray spectrum de- 
tected by Chandra of the ULX X-1 in M82, within 
the framework of Advection Dominated Accretion Flows 
(ADAF) which successfully explains the hard state spec- 
tra of X-ray binaries. They found that if the source is 
likened to the low-luminosity hard state then the black 
hole mass should be M ^ lO^M© else it should be 
~ IO^Mq if the system is to be compared with the high 
luminosity hard state. This degeneracy occurs because 
the observed spectrum during an off axis observation by 
Chandra ( the source is affected by count pile-up for 
normal observations) was a featureless power-law as it 
should be if it is analogues to the hard state. XMM ob- 
servations of the so urce reveal a more coni plex turnover 
at around ~ 8 keV (jAgrawal fc Misral ^2006l , which could 
either be because the spectral state was different during 
the XMM observation or that there were serious contam- 
ination from nearby sources. However, IStobbart et al.l 
(2006) report that many ULX observed by XMM-Newton 
do show such high energy turnovers, which makes the 
comparison with the hard state of black hole binaries 
ambiguous. 



Senorita Devi et al. 



A more robust argument would be provided if a ULX 
reveals soft state like spectral property. In the pure soft 
state of X-ray binaries, the luminosity is dominated by 
emission from a standard disk that extends to the last 
stable orbit. During this state, the high energy power- 
law component (which is dominant during the low and 
intermediate states) contributes < 5% of the total lu- 
minoisty. The state occurs when the Eddington ratio 
> 0.02. During luminosity variations (typically for X- 
ray novae) the constancy of the inner most radius and 
it's value being always close to the predicted last sta- 
ble orbit is taken as strong evidence for the correct- 
ness of the model. In fact, based on the model, re- 
cent attempts have been made to estimate the spin of 
black holes by measuring the subtle strong gravity ef- 
fects of light bend ing and red-shift on the spectra (e.g. 
iMcClintock et all [2006). The analogy to ULX predicts 
that they should exhibit a similar disk black body com- 
ponent with a smaller inner disk t emperature, due to 
the larger black hole mass. Indeed, iMiller et al\ ([20031 ) 
report the presence of a power-law spectrum with a cool 
accretion disk component {kTin ^ 0.1 — 0.5 keV) in NGC 
1313 X-1 and X-2 corresponding to black hole masses 
M ~ W^Mq. a comparison of ULX which show such 
soft components wi th X-ray binarie s has been undertaken 
by [Miller, Fabian fc Milleii (|2004[ ) who argue that the 
systematic lower temperatures indicate that the systems 
harbor IMBH, but also caution for potential weakness of 
such interpretations. The primary issue is the possibil- 
ity of incorrect spectral modeling especially of data with 
low counts ([Goncalvcs & Soria 2006). The soft compo- 
nent is strongly effected by absorption and errors in the 
estimation of the column density may crucially affect the 
results. Another perhaps inconsistent aspect is that for 
these systems, the power-law component is nearly as lu- 
minous as the disk one and hence an analogy has to be 
made with the rarer Very High state (VHS) rather than 
the more frequent pure soft state. Moreover, since the 
power-law component is an important contributer to the 
flux, the modeling of the disk component is suspect to 
the uncertainties in modeling the power-law one. 

The existence of ULX in different spectral states may 
also be revealed in a systematic spectral study of a 
large sample of po tential sources ([Winter et al.l l2006f ) . 
ISwartz et al\ ([2004D fitted the spectra of such a large 
sample with an absorbed power-law model and found no 
bimodal distribution on the spectral index. However, us- 
ing a disk black body model, a bimodal distribution was 
revealed at least for high lumin osity sources for sam ples 
obtained from XMM-New ton ([Winter et al.l l2006f ) and 
Chandra ([Devi et al.ll2007f ). One set of high luminosity 
sources have temperatures kT ^0.1 keV while the other 
have systematically higher temperatures kT ~ 1 keV. 
While the higher temperature sources do not seem to 
have an analogy with black hole binaries, the lower tem- 
perature one may be identified as the equivalent of the 
soft state. For a sample consisting of sources observed by 
Chandra which were not affected by pile-up and which 
were not located in regions of excessive diffuse emission, 
[Devi et aLl ([20071 ) fitted both an absorbed power-law and 
a disk black body spectral model to ascertain the depen- 
dency of luminosity on the spectral model used. They 
identified a highly luminous source, X-7, in NGC 6946 as 
having a soft spectrum and with sufficiently high counts 





TABLE 1 
Spectral properties of NGC 6946, X-7 






kTir, (keV) 


/?i„(xlO«cm) L(x 10^9 erg/sec) 


xVd.o.f 


HCR 
LCR 


29+°"3 

9fi+003 
'J-^'^-0.03 


e.ollj 3.7lg:i 

6.0ll:t 2.81^4 


10.9/16 
22.0/14 



Note. — Best fit spectral parameters for the high rate segment 
(HCR) and the low rate segment (LCR). The model is an absorbed 
disk black body emission with the column density fixed at the Galactic 
value 2.1 X 10 cm^ . The inner radius of the disk B.in is computed 
from the normalization of the disk black body component using the 
distance to the source D — 5.5 Mpc, the viewing angle cosi — 0.5 and 
color factor / — 1.7. 



for more detailed studies 

In the next section we report on spectral and temporal 
analysis of the Chandra observation of this source and 
summarize and discuss the results in the last section. 

2. SPECTRAL AND TEMPORAL PROPERTIES OF NGC 
6946, X-7 

While the temporal beh avior of all the sources exam- 
ined by [Devi et aU ([20071 ) will be presented elsewhere, 
here we report on the ksec variability observed for source, 
X-7, (R.A: 20 35 0.13 and Dec: -1-60 9 7.97 ) in NGC 6946. 

dLh ■ 



The s ource is a known variable source ([Liu fc Mirabell 
[20051: [Lira et al.| 120001) and has been calle d IXO 85 
f Colb ert fc PtakI |2002[) and source no. 56 ([Holt et all 
2003). 

NGC 6946 was observed by Chandra for an effective 
time of 59 ksec on 7 September, 2001 (Observation ID: 
1043). The data reduction and analysis were done using 
CIA03.2 and HEASOFT6.0.2. Using a combination of 
CIAO tools and calibration data, the source (and back- 
ground) spectrum and light curve were extracted. Spec- 
tral fitting was done over the energy range 0.3 to 8 keV 
and spectra were rebinned such that each bin had a min- 
imum of 40 counts. The distance to the source has been 
estimated to b e between 5.1 Mpc |de Va ucoulcurs 197!|) 
and 5.9 Mpc ([Karachentsev et al.l l200(i( ) and hence we 
adopt here a distance of 5.5 Mpc. 

The light curve binned over 2 ksec is shown in Figure 
1, where a clear transition from high to low count rate is 
visible. The probability that the count rate was a con- 
stant during the observations is < 2 x 10"^". The bottom 
panel of the Figure shows the hardness ratio versus time 
binned over 8000 sees. The ratio is defined as the ratio 
between flux in the 0.3 — 1.0 keV band over that in the 
1.0 — 8.0 kcV band. There is evidence that the ratio de- 
creases with the intensity. The two vertical lines mark 
the time ranges used for spectral analysis, with the first 
22.46 ksec being called high count rate segment (HCR) 
state and the last 29.54 ksec as the low count rate one 
(LCR). The total number of counts for the HCR segment 
was 820 while for the LCR segment it was 725. 

The Galactic hydrogen column density along the di- 
rection of the source is 2.1 x 10^^ cm^^. An absorbed 
power-law model with this value of absorption gives an 
unacceptable fit to both HCR and LCR segments, with 
reduced xt > 2.5. The fit can be significantly improved 
by allowing for intrinsic absorption, Nh ^ 10^^ cm^^ but 
with a large unphysical photon index of F > 5. On the 
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Fig. 1. — Top panel: Light Curve of source X-7, binned over 
2000 sees. A clear decrease of counts by a factor ~ 1.5 is visible. 
The vertical lines mark the time which is used to divide the light 
curve into high count rate (HCR) and low count rate (LCR) states 
for spectral analysis. Bottom panel: Hardness ratio versus time 
binned over 8000 sees. The ratio is defined as the ratio between 
flux in the 0.3 — 1.0 keV band over that in the 1.0 — 8.0 keV band. 
There is evidence that the ratio decreases with the intensity. 
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Fig. 2. — The fitted unfolded spectra for the high count rate 
segment, HCR, for parameters listed in Table 1. The solid line 
is the best fit model, while the dashed line represents unabsorbed 
spectra. 



Fig. 3. — The fitted unfolded spectra for the low count rate 
segment, LCR, for parameters listed in Table 1. The solid line is 
the best fit model, while the dashed line represents unabsorbed 
spectra. 



other hand, as shown in Table 1 and illustrated in Fig- 
ures 2 and 3, an absorbed disk black body model with 
the column density at the Galactic value, provides ad- 
equate fits to both segments. The inner radius of the 
disk Rin is computed from the normalization of the disk 
black body component using the distance to the source 
D = 5.5 Mpc, the viewing angle cos i = 0.5 and color 
factor / — 1.7. The result and it implications are better 
represented in Figure 4, where confidence region ellipses 
are shown for the two parameters Rm and inner disk tem- 
perature Tin- The Figure reveals that the simplest and 
more favored interpretation is that the inner disk radius 
remains constant, while the temperature decreases dur- 
ing the variation. Assuming that Rin ~ lOGM/c^, the 
black hole mass, M, can be estimated and confidence el- 
lipses for M and total bolometric luminosity are shown 
in Figure 5. A black hole mass of M ^^ 400Mq is favored 
for both sets while the luminosity decreases from ~ 3.7 
to ^ 2.8 X 10'^^ ergs/sec, corresponding to an Eddington 
ratio, L/LEdd ^ 0.06. 

It is worth emphasizing that the above interpretation 
is the simplest one which provides acceptable fit to the 
data and requires the smallest black hole mass. If one 
allows for intrinsic absorbing column density, the best 
fit temperature decreases and the inner radius increases 
implying a larger black hole mass. For example, a sim- 
ple model could be that the variation is only due to 
changes in the intrinsic column density. A model where 
the column density is allowed to vary, with Rin and 
Tin same for both segments, indeed gives an acceptable 
combined x^/dof = 36.9/32, with Nh increasing from 
3 ± 0.4 to 3.9 ± 0.4 X 10^^ cm"^ and inner disk tem- 
perature Tin =■ 0.22 ±0.1 keV. The inner disk radius 
turns out to be 1.5 ± 0.3 x lO'' cm corresponding to a 
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Fig. 4. — Confidence region eUipses for the inner disk radius, Ri„ 
and temperature, Ti„, for the HCR and LCR segments for best fit 
parameters given in Table 1. The soHd (dashed) curves represent 
Ax^ = 2.3(4.61) corresponding to 68.3 (90%) confidence level. The 
results favor the interpretation that Rin is a constant and that the 
decrease in count rate is due to a decrease in Ti„. 
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Fig. 5. — Confidence region ellipses for the black hole mass, 
M and bolometric luminosity, L, for the HCR and LCR segments 
estimated from iJi„ and Ti„ values shown in Figure 3. The solid 
(dashed) curves represent Ax^ = 2.3(4.61) corresponding to 68.3 
(90%) confidence levels. For both the segments a black hole mass, 
M ~ 4OOM0 is favored. 



~ 95OAf0 black hole. If the column density and tem- 
perature are allowed to vary independently for the two 
segments, x^/dof = 29/29, with i?„ = 1.3 ± 0.4 x 10^ 
cm corresponding to M ~ 900Mq. An additional power- 
law component with photon index T = 2.0 or F = 2.7 
does not improve the fit if the column density is fixed 
at the Galactic value. For Nh '^ 6 x 10^^, a disk black 
body and a power-law component with F = 2, provides 
a good fit to both segments x^/dof = 19.9/27. Here, 
the variability is due to a change of normalization of the 
power-law component, while the temperature Tm ~ 0.13 
keV and i?i„ ~ 1.2 x 10^" cm remain nearly constant. 
Again, in these more complicated interpretations Rin is 
significantly larger implying a larger black hole mass. 

3. SUMMARY AND DISCUSSION 

The X-ray source X-7 in NGC 6946 is located in a 
region of low diffuse emission and has a count rate which 
does not cause pile-up in Chandra detectors. The light 
curve of this source reveals a decrease of ~ 1.5 in the 
count rate over 5000 sec making it one of the few ULX 
that have clearly shown variability on ksec time-scales 
(jKrauss et al.ll2005t iMizuno et al.ll2001h . 

Spectral fitting of the high and low count rate parts of 
the light curve reveals that a simple model of a disk black 
body emission absorbed by the Galactic column density 
of Nh — 2.1 X 10^^ cm~^ can adequately represent both 
segments. The best fit temperature varies from 0.26 to 
0.29 keV while the inner disk radius remains constant at 
'--^ 6 X 10* cm. This would imply that the mass of the 
black hole is r^ 400Mo and at a luminosity of ~ 3 x 10'^^ 
ergs/sec the system has an Eddington ratio of ~ 0.06. 
Other more complicated spectral fits like assuming in- 
trinsic column density for the source and/or addition of 
a power-law component results in a larger estimation of 
inner disk radius and hence a larger black hole mass. 
Although the luminosity of the source estimated here 
is only mildly super-Eddington for ten solar mass black 
hole, the low inner disk temperature ~ 0.3 keV implies a 
larger inner disk radius and hence a large black hole mass 
~ 400Mo. A low sub-Eddington accretion rate on a ten 
solar mass black hole could produce the observed low in- 
ner disk temperature, but then the predicted luminosity 
would be sig nificantly less than wha t is ob served. As 
cautioned bv iMiller. JFabian fc Milled ()2004l ). the source 
could be a background AGN and the soft component ob- 
served is actually the soft excess which is detected in 
many AGN. These soft excesses can be modeled as black 
body emission at kT ~ 0.1 keV similar to the soft com- 
ponent in ULX. While this may be true for some ULX, 
this is unlikely in this case because here (unlike soft ex- 
cess of AGN) the soft component totally dominates the 
luminosity. 

If the emission is beamed (jKing et aLll2001h by a factor 
7], the inner disk radius would be overes timated by 77^/^. 
Howev er, even for extreme geometries iMisra fc SrirairJ 
()2003f ) have computed rj < 5. Hence an extreme beam- 
ing of 77 ~ 5, would imply that Rm "^ 2.7 x 10* cm 
corresponding to a black hole mass of ISOMq. The color 
factor used in this analysis is the standard value of 1.7. 
Even for the extreme case that there is no color correc- 
tion, Rin ^ 2x 10* cm corresponding to a black hole mass 
of 130AfQ. Thus only in the extremely unlikely case of 
having no color correction and high beaming factor 77 ~ 5 
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would the black hole mass be < IOOMq. 

A model can be envisioned where the inner disk radius 
is not at lOGM/c^ but is at a much larger radius say 
lOOGM/c^, with the inner region being highly radiatively 
inefficient. Then the mass of the black hole could be 
~ 40Mq. However to produce a luminosity of ^^ 2 x 10'^^ 
ergs/sec at lOOGM/c^ would require an accretion rate 
which is 50 times the Eddington value and it is not clear 
how the inner region would remain radiatively inefficient 
at such accretion rates. A more radical model would be 
that the emission arises from an optically thick sphere of 
radius ^ lOOGM/c^ surrounding a AQMq black hole, but 
that would be a serious paradigm shift from our present 
understanding of accretion flows. 

In a contemporary work, Fridri ksson et a l. (2008"), 
have studied the long term variability of sources in NGC 
6946 using five Chandra observations. They confirm the 



variability of this source (their source no. 60) for this 
observation, but not for the others. They detect long 
term flux variability for the different observations. Their 
hardness ratio estimations confirm the soft nature of this 
source for all the observations. 

The results presented here indicate that the black hole 
mass of this source is > 400Mq which may be true for 
other super-soft X-ray sources and for ULX in general. 
However, further analysis of data of this and other super- 
soft sources are required to ascertain whether such tem- 
poral and spectral behavior is common among ULX. 
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